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\ ^ • Abstract 

We present perfect fluid Friedmann-Robertson- Walker quantum cosmological models in the presence of 
negative cosmological constant. In this work the Schutz's variational formalism is applied for radiation, 
dust, cosmic string, and domain wall dominated Universes with positive, negative, and zero constant spatial 
curvature. In this approach the notion of time can be recovered. These give rise to Wheeler- De Witt 
equations for the scale factor. We find their eigenvalues and eigenfunctions by using Spectral Method. 
After that, we use the eigenfunctions in order to construct wave packets for each case and evaluate the 
time-dependent expectation value of the scale factors, which are found to oscillate between finite maximum 
6JT), and minimum values. Since the expectation values of the scale factors never tends to the singular point, 

we have an initial indication that these models may not have singularities at the quantum level. 
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1 Introduction 



Quantum cosmological models are important subjects on the interface of cosmology and gravitation. At first, 
B. DeWitt [I] quantized a Friedmann Universe filled with dust and later, closed isotropic cosmological models 



with matter as a conformal and minimally coupled scalar fields were quantized [2l [3]. Misner worked on 
the quantization of anisotropic cosmological models [3], and Barabanenkov quantized the Friedmann metric 
matched with the Kruskal one The quantization of a dust-like closed isotropic cosmological model with a 
cosmological constant is also investigated in Ref. [BJ. 

In the quantum cosmology the Wheeler-DeWitt (WD) equation which determines the wave function of the 
Universe, can be constructed using ADM decomposition of the geometry [7] in the Hamiltonian formalism of 
general relativity. However, quantum cosmology has many technical and conceptual problems. In fact, the 
WD equation of quantum gravity is a functional differential equation defined in the superspace which is the 
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space of all possible three dimensional spatial metrics, and no general solution is known in this superspace. In 
quantum cosmology this problem is avoided by using symmetry requirements to freeze out an infinite number 
of degrees of freedom, leaving only a few for quantization process. This procedure defines a minisuperspace, 
where exact solutions can often be found. On the other hand, the general covariance will be lost upon applying 
the ADM decomposition, and in most cases the notion of time disappears at the quantum level [5j. Even, if 
all these problems are solved, the interpretation of the central object, i.e. the wave function of the Universe, 
remains unanswered. 

The many- worlds interpretations [9 J of quantum mechanics is one of the most popular interpretation schemes 
for the wave function of the Universe. This interpretation differs noticeably from the Copenhagen interpretation 
of quantum mechanics since the conception of probability is abandoned in some sense. In fact, all possibilities 
are participated to create new Universes with different possible eigenvalues obtained by measurements. The 
evolution of observables such as scale factor is found by evaluating the expectation values. In this case, like in 
the Copenhagen interpretation, the structure of Hilbert space and self-adjoint operators are still unchanged. 

The presence of the matter in quantum cosmology needs further consideration and can be described by 
fundamental fields, as done in Ref. [10]. Using WKB approximation one can predict the behavior of the 
quantum Universe which leads to determination of the trajectories in phase space. However, even in the 
minisuperspace, general exact solutions are hard to find, the Hilbert space structure is ambiguous and it is 
difficult to recover the conception of a semiclassical time [U [10] . 

Here, we consider matter as a perfect fluid. This description is basically semiclassical, but it introduces a 
variable, which can be identified with time and connected with the matter degrees of freedom, leading to a 
well-defined Hilbert space structure. Moreover, this allow us to treat the barotropic equation of state p = ap 
with arbitrary a. 

It is very convenient to construct a quantum perfect fluid model. Schutz's formalism 1X2] gives dy- 
namics to the fluid degrees of freedom in interaction with the gravitational field. Using a proper canonical 
transformations, at least one conjugate momentum operator associated with matter appears linearly in the 
action integral. Therefore, a Schrodingcr-like equation can be obtained where the matter plays the role of 
time. Moreover, recently, some applications of the Schutz's formalism have been discussed in the framework of 
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the perfect fluid Stephani Universe [T3l [14] and Friedmann-Robertson- Walker (FRW) Universe in the presence 
of Chaplygin gas [HI [16] . 

Until now, quantum perfect fluid models with common equations of state have been constructed in the 
absence of cosmological constants jTTJ [TBJ [TU [501 HI] • We can study the behavior of the scale factor using the 
many-worlds and the de Broglie-Bohm interpretations of quantum mechanics. 

Recently, the quantization of FRW radiation dominated Universe in the presence of a negative cosmological 
constant is discussed by Monerat et al |22| . However, as mentioned in Ref. [23], their results are inaccurate and 
their relative errors range between 10 -3 for the ground state of k = 1 case, and 1 for the ground state of k = — 1, 
which make their work unreliable. Here, we generalize the previous investigations by studying quantum perfect 
fluid models for barotropic equation of state p = ap, where a = {1/3,0, —1/3, —2/3} correspond to radiation, 
dust, cosmic string, and domain wall dominated Universes, respectively. Using the many-worlds framework, 
the behavior of the scale factor is determined, although the results are independent of the interpretation scheme 
employed. The large time average of the expectation value of the scale factor is similar to the classical case. 
Moreover, the model predicts an accelerated expansion today if— l/3>a>— 1. 

It is important to mention that although recent observations point toward a positive cosmological constant, 
it is still possible that at the very early Universe the cosmological constant be negative. Moreover, we think it 
is important to understand more about such models which represent bound Universes. 

This paper is organized as follows. We quantize three Friedmann-Robertson- Walker perfect fluid models in 
the presence of a negative cosmological constant, using the formalism of quantum cosmology. In Sec. [U the 
quantum cosmological model with a perfect fluid as the matter content is constructed in Schutz's formalism 

I12j. and the WD equation in minisuperspace is Found to quantize the model. The wave-function depends 
on the scale factor a and on the canonical variable associated to the fluid which plays the role of time T, in the 
Schutz's variational formalism. We separate the wave-function in two parts, one depending solely on the scale 
factor and the other depending only on the time. The solution in the time sector of the WD equation is trivial, 
leading to imaginary exponentials of the type e , where E is the system energy and t = T. In Sec. [3] we 
outline the Spectral Method pBJ [2SJ [26] , and use it to find the eigenvalues and eigenfunctions of corresponding 
WD equations. In Sec. 01 we construct wave packets from the eigenfunctions, for radiation, dust, cosmic string 
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and domain wall dominated Universes respectively, and compute the time-dependent expectation values of the 
scale factors for k = 1,0,— 1. In Sec. El we state our conclusions. 

2 Model 

Let us start from the Einstein- Hilbert action plus a perfect fluid in the formalism developed by Schutz. For 
this, we write down the action for gravity plus perfect fluid as 

S = - I d 4 x^{R- 2A) + 2 f d 3 xVhh ab K ab + [ d 4 x</=gp, (1) 

2 Jm JdM JM 

here, K ab is the extrinsic curvature, A is the cosmological constant, and h a j, is the induced metric over the 
three-dimensional spatial hypersurface, which is the boundary DM of the four dimensional manifold M. We 
choose units such that the factor 8ttG becomes equal to one. The first two terms were first obtained in [7] and 
the last term of (|TJ) represents the matter contribution to the total action, p being the pressure which obeys the 
barotropic equation of state p — ap. In Schutz's formalism |11[ 1 1 2 j the fluid's four-velocity can be expressed 
in terms of five potentials e, £, /3, 9 and S 

u„ = -{£,» + (P,v + 0S,u) (2) 
M 

where p, is the specific enthalpy. S is the specific entropy, and the potentials £ and (3 are connected with 
rotation which are absent of models in the FRW type. The variables e and 9 have no clear physical meaning. 
The four-velocity also satisfies the normalization condition 

u v u u = — 1. (3) 

The FRW metric 

ds 2 = -N 2 (t)dt 2 + a 2 (t)g l3 dx l dx\ (4) 

can be inserted in the action |T]), where N(t) is the lapse function and gij is the metric on the constant- 
curvature spatial section. After some thermodynamical considerations and using the constraints for the fluid, 
and dropping the surface terms, the final reduced action takes the form [18 . 



S 
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The reduced action may be further simplified using canonical methods [18] , resulting in the super-Hamiltonian 



H 



12a 



Aa 3 - 3ka + 



n 3a ' 



(6) 



where p a = — 6aa/7V and p e — —pou°Na 3 , po being the rest mass density of the fluid. The following additional 
canonical transformations, which generalizes the one used in [18] , 



T = 



- Ps e~ S p7 ia+1) , 



I | -i \ P s 
e- (a + 1 — , 

Pe 



PT = p« +1 e s 



Pe =Pe, 



(7) 



simplifies the super-Hamiltonian to, 



n 



Pi, , a 3 o; , PT 

+ Aa — Ska + 



12a 



3n 



(8) 



where the momentum px is the only remaining canonical variable associated with matter and appears linearly 
in the super-Hamiltonian. The parameter k defines the curvature of the spatial section, taking the values 
0, 1, —1 for a flat, close or open Universes, respectively. 

The classical dynamics is governed by the Hamilton equations, derived from Eq. ([5]) and Poisson brackets 



as 



Np a 
6 a 



{ Pa ,NH} = -j^pl + 3Nk 



(9) 



a= {a 7 Nn} 
Pa ■ 

-37V Aa 2 + NZaar Za - 1 pT , 
f = {T, NH} = Na- 3a , 
Pt = {pt, NH} = . 

We also have the constraint equation 7i = 0. Choosing the gauge N = ait), we have the following solutions 
for the system 



-ka + -Aa 3 
3 



1 — 3a 



G 



a' 3a p T , 



= -3a 2 + Aa 4 - 3ka 2 + a 1 ' 3 " 



Pt- 



(10) 
(11) 



The classical equation of motion for the scale factor in absent of the cosmological constant is solved in a unified 
form for any a € [0, 1] in terms of hypergeometric functions in Ref. |27j . Moreover, In the radiation dominated 



Universe (a = 1/3) with a negative cosmological constants, the classical solutions have been obtained using 
Jacobi's elliptic sine functions |22j . the WD equation in minisuperspace can be obtained by imposing the 
standard quantization conditions on the canonical momenta and (p a — > — Pt — ► — demanding that 
the super-Hamiltonian operator annihilate the wave function (Ti = 1) 

+ 12Aa% - 36ka 2 »£ - iUa 1 ' 301 — = 0. (12) 
oa z at 

where t — T corresponds to the time coordinate. Equation (|12p takes the form of a Schrodinger equation 
id"i>/dt = H^. Demanding that the Hamiltonian operator H to be self-adjoint, the inner product of any two 
wave functions $ and \& must take the form [28l [19] 

/>oo 

($,*) = / a 1 - 3a $*1'cfc, (13) 
On the other hand, the wave functions should satisfy the following boundary conditions 



*(0,<) = or 



= 0. (14) 

a=0 



da 

The WD equation (fT2")) can be solved by separation of variables as follows 

i/j(a,t)=e im i>(a), (15) 

where the a dependent part of the wave function (ip(a)) satisfies 

- ip"(a) + (36fca 2 - 12Aa 4 )V'(a) = 12Ea 1 ' 3a ij(a). (16) 

Since the energy term grows faster than the potential for a < — 1, this equation has a discrete spectra (E n ) 
with associated bound state eigenfunctions (ip n (x)) only for a > — 1. 

We construct a general solution to the WD equation (fT2")) by taking linear combinations of the - 0n( a ! ^)' s > 

m 

*(o, t) = J2 C n (E n )iP n (a)e tEnt , (17) 

n=0 

where the coefficients C n {E n ) will be fixed later. Form pure mathematical point of view, by allowing negative 
values of a, the Parity operator can be defined. If the WD equation (Eq. [T5|) is covariant under the Parity 
operator, its eigenfunctions can be separated into even and odd ones. The even or odd wave packets constructed 



from appropriate linear combinations of the eigenstates, have the important property that they will not change 
their parity in the course of their time evolution. Therefore, if we choose the initial wave packets to be odd or 
even, that is they satisfy either the first or the second condition stated in Eq. [HI respectively, they will satisfy 
them for all times. We can compute the expected value for the scale factor a for any wave function, using the 
many worlds interpretation of quantum mechanics. This means, we may write the expected value for the scale 
factor a as [29] 

_ io°° a 2 - 3a \^( ai t)\ 2 da 
{a)t ~ J™a^\*(a,t)\ 2 da (18) 

Before solving the WD equation (fTr?)) via Spectral Method, it is worthy to state a brief overview of the 
Chhajlany and Malnev method and Variational Sine Collocation Method (VSCM), which have been recently 
used to solve the WD equation (fT5]) . for radiation epoch (a — 1/3) in Refs. [221123], respectively. 

In Chhajlany and Malnev method [301 OS], one adds an extra term to the original anharmonic oscillator 
potential to find a subset of normalizable solutions of the modified Hamiltonian. In the case of equation (|16| 
this extra term is ca e where c is constant. Now, the solution can be written as a polynomial where the larger 
the degree of the polynomial, the smaller the constant, c is. In fact, by increasing the order of polynomial, the 
energy eigenvalues predicted by this method approach monotonically to the energy eigenvalues of the original 
Hamiltonian. 

On the other hand, to obtain highly accurate numerical results, both for the energy eigenvalues and eigen- 
functions, one can use Variational Sine Collocation Method (VSCM) [32]. It is shown that the errors decay 
exponentially with the number of elements (sine functions) used for discretization of the Hamiltonian. Diago- 
nalization of the resulting matrix, by specification of the otherwise arbitrary grid spacing h (spacing between 
two contiguous sine functions), yields energy eigenvalues and eigenfunctions. As shown by Amore et al, for a 
specified number of sine functions, there exists an optimal value of h which yields the minimum errors [32] , 
this optimal value can be found using the Principle of Minimal Sensitivity (PMS) [33 to the trace of the 
Hamiltonian matrix. 

As indicated by Lemos et al, the need for a modified potential instead of the original one in the Chhajlany 
and Malnev method, gives rise to significant errors, particularly for k = — 1 [34]. In fact, VSCM is more 
uniformly accurate and converges more rapidly than the Chhajlany and Malnev method. 



3 The Spectral Method 

In this section we introduce Spectral Method (SM) J3H[25] as a tool for solving differential equation. We have 
recently used the Spectral Method for constructing the appropriate wave packets which are solutions to a WD 
equation [26]. This method is simple, fast, accurate and stable. 

Let us consider the general time- independent WD equation (EQ. (fl6| ). 



dx 2 



+ f[xW(x)=Eg[x]Mx), (19) 



where / and g are arbitrary, but with derivative operators less than two. For the usual eigenvalue problem 
g = 1, which includes the time-independent Schrodinger equation. The method SM can be easily extended 
to solve the general case which g is a operator in the x space. This generalize problem can be named a 
generalized eigenvalue problem. Throughout this paper, we only examine the bound states of this problem, 
i.e. the states which are the square integrable. The configuration space for most physical problems are defined 
by — oo < x < oo. Since the bound states fall off sufficiently fast for large \x\, a finite region suffices, and 
the proper choice for this region, say — L/2 < x < L/2. The use of a finite domain is also necessary since we 
need to choose a finite subspace of a countably infinite basis. We find it convenient to shift the domain to 
< x < L. In particular, we need to shift the potential energy functions also. This means that we can expand 
the solution as, 



</>(x) = 5> n J- sin (— )• (20) 

71=1 

We can also make the following expansions, 

U(x) = y|sin(^), (21) 

n 

g i){x) = J2 B n y 1 « n (— ) . ( 22 ) 



where B n B' n are coefficients that can be determined once / and g are specified. By substituting Eqs. (|20I21I22[) 
into Eq. (|19p and using the differential equation of the Fourier basis we obtain, 

E/n7r\ 2 . „ 1 . fnirx\ „, . /nirx\ . „. 

(—J A n +B n sin(— J =EJ2K^{— )■ (23) 



Because of the linear independence of sin (^r^j, every term in the summation must satisfy, 

Y A n + B n = E B' n . (24) 



/ T17T ^ 



\ L J 

It only remains to determine the matrices B and B 1 . Using Eqs. (|21I22|) and Eq. (|2"0]) we have, 

n n 

K sin (— ) = £ A n £ sin (_) , (26) 

n n 

By multiplying both sides of the above equations by sin and integrating over the x-space and using the 

orthonormality condition of the basis functions, one finds, 

m 

B 'n = X! ^.n ^-to, (28) 



where, 



2 ' L 



r' = 



/ sin(=)/ S in(^)^ (29) 
(^)^sin(^)^ (30) 



I / . rrmx 



Therefore we can rewrite Eq. ([24]) as, 

+ ^ ] Crn.n A m = E ^ ' C m n A m . (31) 

It is obvious that the presence of the operators / and g in Eq. (TIT))) , leads to nonzero coefficients C m .„ and 
C' m n in Eq. (f3"Tj) , which in principle could couple all of the vector elements of A. It is easy to see that the 
more basis functions we include, the closer our solution will be to the exact one. We select a finite subset of 
the basis functions i.e. the first N ones, by letting the index m run from 1 to N in the summations. Equation 
(f3Tj) can be written as, 

DA = ED'A, (32) 



or. 



D'- 1 DA = EA, (33) 
9 



where D and D' are square matrices with N x N elements. Their elements can be obtained from Eq. (|3ip . 
The solution to this matrix equation simultaneously yields N sought after eigenstates and eigenvalues. It is 
important to note that the optimized value of L crucially depends on the number of basis functions N (L(N)), 
which results in the maximum accuracy and the stability of the solutions (for a comprehensive study about 
the optimization procedure see [25 ). 

4 Results 

In this section we will solve the equation ([T5)l using SM. By choosing N = 100 basis functions, and we report 
our results with 10 significant digits. Note that, although, we are free to choose other values of A, but the 
accuracy of results for small |A| reduces in comparison with large values of |A| for a given number of basis 
TV, particularly for k = — 1. This means that we need to increase the number of basis N to obtain the same 
accuracy which increases the computations. With regard to these considerations, the results are robust under 
changes of Lambda. 

4.1 Radiation (a = 1/3) 

In the radiation dominated Universe time-independent WD equation has the following form, 

d 2 ip{a) 



da 2 



+ (36ka A - 12Aa 4 )?/>(a) = 12Eip(a), (34) 



In this form it is obvious that the system is absolutely stable for A < 0. Note that equation is covariant under 
the Parity operator. For ease of comparison of our results with those of Refs. [221 [23] , we select the first 
condition of the equation (|14|) and choose the coefficients C„ s in equation (|17|) to be 1 and zero for the odd 
and even eigcnfunctions, respectively. We can find the energy eigenvalues and eigenvectors of this equation 
with ease using SM where / = 36kx 2 — 12Aie 4 and g — 12 in comparison with Eq. ([HI). Table Q] shows the 
first 26 odd eigenvalues for k = 1,0,-1 respectively. Figures i] 1 12l;j|) show the resulting expectation values of 
the scale factor a, versus t for the various values of k. As can be seen from the table, the results are as same 
as those reported in Ref. [23j . To show the arbitrariness in choosing initial odd wave packet, we can use the 
coefficients of odd coherent state of the quantum simple harmonic oscillator. Figure [4] shows the 3D plot of 
resulting wave packet for k = 1 case. 
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Figure 1: Behavior of the expectation value of the scalar factor for A = —0.1, k = 1, and C n = 1,0 for odd 
and even n, respectively, in radiation regime. 
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Figure 2: Behavior of the expectation value of the scalar factor for A = —0.1, k = 0, and C n = 1,0 for odd 
and even n, respectively, in radiation regime. 
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Figure 3: Behavior of the expectation value of the scalar factor for A = —0.1, k = —1, and C n = 1,0 for odd 
and even n, respectively, in radiation regime. 
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Figure 4: 3D plot of the square of the wave packet for fc = 1 in radiation regime with the coefficients of odd 
coherent state of the quantum simple harmonic oscillator. 
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E x 


1.510262538 


0.3364795921 


-21.79569604 


E 2 


3.550647291 


1.031199050 


-20.39848969 


E 3 


5.621893706 


1.880761581 


-19.01885126 


E 4 


7.722777814 


2.842487493 


-17.65745105 


E 5 


9.852194220 


3.894746211 


-16.31503013 


E 6 


12.00913857 


5.024091561 


-14.99241327 


E 7 


14.19269339 


6.221192430 


-13.69052532 


E s 


16.40201655 


7.479120856 


-12.41041248 


E 9 


18.63633173 


8.792490724 


-11.15327022 


Eio 


20.89492047 


10.15697184 


-9.920481041 


En 


23.17711552 


11.56899282 


-8.713666752 


E\2 


25.48229516 


13.02554815 


-7.534763378 


E13 


27.80987836 


14.52406697 


-6.386132854 


Eu 


30.15932069 


16.06232038 


-5.270738293 


E15 


32.53011067 


17.63835396 


-4.192437789 


Ei6 


34.92176669 


19.25043737 


-3.156519423 


En 


37.33383429 


20.89702582 


-2.170719225 
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39.76588373 


22.57673016 


-1.246288323 


E19 


42.21750793 


24.28829319 


-0.3899963301 
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44.68832056 


26.03057067 


0.4337198672 


E21 


47.17795441 


27.80251593 


1.300741394 


E22 


49.68605987 


29.60316715 


2.243855391 


E23 


52.21230364 


31.43163692 


3.256001774 


E 2 4 


54.75636750 


33.28710339 


4.326509414 


E 25 


57.31794728 


35.16880284 


5.448058248 


E26 


59.89675181 


37.07602341 


6.615611361 



Table 1: The lowest calculated energy levels for the cases fc = 0, fc = 1, and k = — 1 in radiation dominated 
Universe (in all cases, A = —0.1). 
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4.2 Dust (a = 0) 

In dust dominated Universe time-independent WD equation has the following form, 

- + (36ka 2 - 12Aa 4 )^(a) = UEaMa), (35) 

da 2 

We can find the energy eigenvalues and eigenvectors of this equation with ease using SM where / = 36kx 2 — 

12Ax 4 and g — 12x in notation displayed in Eq. (|19p. Table [2] shows the first 20 positive eigenvalues for 

k = 1,0,-1 respectively. Note that, for any positive eigenvalues (E+), there is an negative counterpart (E~) 

which E~ — —E+. The above equation, is not invariant under the Parity operator. Therefore, its eigenfunctions 

can not in general satisfy either of the conditions of equation (|141 However, we can construct wave packets, 

from linear combinations of the eigenfunctions, which vanishes at a — and t — 0. Then we need to check that 

the constraints (Eq. [TJ] remain valid for all t for our choice of initial condition. For example we can choose 

the coefficients C n s so as to construct a gaussian initial wave packet (^(a,0)) which is centered e.g. at a = 1. 

Figures (|5|6|7p show the resulting expectation values of the scale factor a, versus t for the various values of k. 

As can seen from the figures these wave packets always satisfy the first boundary condition (Eq. IT4|) . 

In the case A = 0, the time- independent WD equation ([33)) reduces to 

_ d 2 j;(a) 
da 2 

In terms of the new variable x = 6a — E we find 



+ 36/ca>(a) = 12Eaip(a), (36) 



d 2 ^P(x) 
dx 2 



+ 



x 2 E 2 



1>(x) = 0, (37) 



36 36 

Equation (|3T[) is formally identical to the time-independent Schrodinger equation for a harmonic oscillator with 
unit mass and energy A: 

- + [-2A + w 2 x 2 ] ((a) = 0, (38) 

where 2A = £' 2 /36 and w = 1/6. Therefore, the allowed values of A are (n + l/2)w and the possible values of 
E are 



E n = ±v/6(2n + l) , n = 0, 1, 2, ... . (39) 

Thus the stationary solutions are 

* n (o, t) = e- <E »Vn (12a - E n ) , (40) 
13 
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Figure 5: Behavior of the expectation value of the scalar factor for A = —15, k = 1, and ^(a,0) = exp(— 8(a- 
l) 2 ) in dust regime. 
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Figure 6: Behavior of the expectation value of the scalar factor for A 
l) 2 ) in dust regime. 



-15, k = 0, and *(a,0) = exp(-8(a- 



where 



<Pn(x) = H n 



-x 2 /24 



(41) 



and H n are the n-th Hermite polynomial. 
4.3 Cosmic Strings (a = —1/3) 

In Cosmic Strings dominated Universe time-independent WD equation has the following form, 

<Pi/>(a) 



da 2 



(36fca 2 - 12Aa 4 )V>(a) = 12Ea 2 ip{a), 



(42) 



This differential equation is covariant under parity operator and hence its eigenstates can be separated into 
even and odd ones. We can find the energy eigenvalues and eigenvectors of this equation with ease using SM 
where / = i&kx 2 — 12Aa; 4 and g — 12x 2 in comparison with Eq. (fT9|) . Table [3] shows the first 20 eigenvalues 
for k = 1, 0, —1 respectively. By choosing the first condition of the equation HH), the resulting wave packets 
should consist of only the odd eigenfunctions. Therefore, the coefficients C n s in equation (fT7|) are arbitrary 
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Figure 7: Behavior of the expectation value of the scalar factor for A = — 15, k = — 1, and $>(a, 0) = exp(— 8(a — 
l) 2 ) in dust regime. 
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53.50529876 


50.31152949 


47.09859109 


Es 


60.35022067 


57.01973343 


53.67089213 


E 9 


67.18479286 


63.72793736 


60.25341729 


Eio 


74.01077406 


70.43614129 


66.84443864 


En 


80.82948903 


77.14434522 


73.44265237 


El2 


87.64196335 


83.85254916 


80.04704775 


E13 


94.44900906 


90.56075309 


86.65682363 


En 


101.2512814 


97.26895702 


93.27133299 


ElB 


108.0493176 


103.9771610 


99.89004490 


Eia 


114.8435643 


110.6853649 


106.5125177 


En 


121.6343972 


117.3935688 


113.1383797 


Eis 


128.4221359 


124.1017728 


119.7673145 


ElQ 


135.2070546 


130.8099767 


126.3990507 


E20 


141.9893905 


137.5181806 


133.0333530 



Table 2: The lowest calculated energy levels for the cases k = 0, k = 1, and k = — 1 in dust dominated Universe 
(in all cases, A = —15). As mentioned in the text for every positive eigenvalue there exist a corresponding 
negative one with identical absolute value. 
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Figure 8: Behavior of the expectation value of the scalar factor for A = — 15, k = 1, and C n = 1,0 for odd and 
even n, respectively, in cosmic strings regime. 




Figure 9: Behavior of the expectation value of the scalar factor for A = — 15, k = 0, and C„ = 1, for odd and 
even n, respectively, in cosmic strings regime. 

for the odd eigenfunctions zero for the even ones. To be able to extend the results of Refs. [22l [23] for the 
radiation case to the present one, we choose the same initial state as their's. That is the odd ones are all chosen 
to be equal to one. Figures (|8I9I10|) show the resulting expectation values of the scale factor a, versus t for the 
various values of k. 

4.4 Domain Walls (a = -2/3) 

In Domain Walls dominated Universe (a = —2/3) the time-independent WD equation has the following form, 



da? 



(36fca 2 - 12Aa 4 )V>(a) = 12Ea 6 ip(a), 



(43) 



We can find the energy eigenvalues and eigenvectors of this equation with ease using SM where / = 36kx 2 — 
12 Arc 4 and g = I2x 3 in comparison with Eq. (I19p . Table 2] shows the first 20 eigenvalues for k = 1,0,-1 
respectively. Note that, for any positive eigenvalues ), there is an negative counterpart (E~) which E~ — 
This case is similar to the Dust case and in particular its differential equation is not covariant under 
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Figure 10: Behavior of the expectation value of the scalar factor for A = —15, k = —1, and C n = 1,0 for odd 
and even n, respectively, in cosmic strings regime. 





k=l 


k = 


k = -1 


Ei 


11.63722935 


8.637229353 


5.637229353 


E 2 


19.36451595 


16.36451595 


13.36451595 


E 3 


25.54276474 


22.54276474 


19.54276474 


E 4 


30.96152618 


27.96152618 


24.96152618 


E 5 


35.89633715 


32.89633715 


29.89633715 


E e 


40.48441592 


37.48441592 


34.48441592 


E 7 


44.80660320 


41.80660320 


38.80660320 


Es 


48.91557892 


45.91557892 


42.91557892 


Eg 


52.84808309 


49.84808309 


46.84808309 


Eio 


56.63102596 


53.63102596 


50.63102596 


En 


60.28486358 


57.28486358 


54.28486358 


El2 


63.82560554 


60.82560554 


57.82560554 


E13 


67.26607937 


64.26607937 


61.26607937 


En 


70.61676389 


67.61676389 


64.61676389 


ElB 


73.88635930 


70.88635930 


67.88635930 


Eia 


77.08218981 


74.08218981 


71.08218981 


En 


80.21050097 


77.21050097 


74.21050097 


Eis 


83.27676241 


80.27676241 


77.27676241 


E w 


86.28668819 


83.28668819 


80.28668819 


E20 


89.25169642 


86.25169642 


83.25169642 



Table 3: The lowest calculated energy levels for the cases k = 0, k = 1, and k = — 1 in cosmic strings dominated 
Universe (in all cases, A = — 15). 
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Figure 11: Behavior of the expectation value of the scalar factor for A = —15, k — 1, and ^(a. 0) = exp(— 8(a — 
1.5) ) in domain walls regime. 

1 . 6\ 
1 . 4 
1.2 

1 

200 400 600 800 1000 

t 

Figure 12: Behavior of the expectation value of the scalar factor for A = —15, k — 0, and ^(a. 0) = exp(— 8(a — 
1.5) 2 ) in domain walls regime. 

Parity Operator and therefore, its eigcnfunctions can not in general satisfy either of the conditions stated in 
equation (|14l However, we can construct wave packets, from linear combinations of the eigenfunctions, which 
vanishes at a = and t — 0. Then we need to check that the constraints (Eq. [14] remain valid for all t for our 
choice of initial condition. For example we can choose the coefficients C„ s so as to construct a gaussian initial 
wave packet (^(a,0)) which is centered e.g. at a — 1.5. Figures (j 11 1 1 21 13[) show the resulting expectation 
values of the scale factor a, versus t for the various values of A:. As can seen from the figures these wave packets 
always satisfy the first boundary condition (Eq. [T4")) . 

It is important to note that we have repeated the simulations for all cases (a = 1/3,0, —1/3, —2/3) with 
other values of A and different initial conditions (subject to ^(0, 0) = 0). In particular, we have also repeated 
simulations for A = —10, —12.5, —17.5, —20 rather than A = —15 which studied in detail, and found the 
corresponding eigenvalues and eigcnfunctions with desired accuracy. Moreover, we chose other initial conditions 
in the form \I/(a, 0) = exp(— 7(0 — ao) s ) with various choices of 7 (2, 5, 10, 20), S (2, 4, 6), and ao (1, 1.2, 1.4, 1.6). 
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Figure 13: Behavior of the expectation value of the scalar factor for A = —15, k = — 1, and *(a, 0) = 
exp(— 8(a — 1.5) 2 ) in domain walls regime. 





k = l 


k = 


k = -1 


Ei 


17.07778092 


12.54649750 


7.458961879 


Ei 


21.41791925 


18.18718412 


14.81583000 


E3 


24.32147857 


21.57780684 


18.75859852 


E 4 


26.60465935 


24.14422756 


21.63352524 


E 5 


28.52359398 


26.25664404 


23.95258205 


Eg 


30.19749314 


28.07460050 


25.92257362 


E 7 


31.69296487 


29.68328712 


27.64978790 


E s 


33.05153398 


31.13420271 


29.19682156 


E 9 


34.30107228 


32.46114131 


30.60396527 


E10 


35.46131513 


33.68761746 


31.89883581 


En 


36.54683096 


34.83073057 


33.10127373 


E12 


37.56911074 


35.90336205 


34.22606529 


E13 


38.54199370 


36.91577767 


35.28463557 


En 


39.50080454 


37.87984925 


36.28735433 


Ei 5 


40.512127 


38.825301 


37.253593 


Eie 


41.627809 


39.812705 


38.232631 


En 


42.859939 


40.895275 


39.288909 


Eis 


44.203790 


42.088418 


40.454402 


•Ei 9 


45.653137 


43.388462 


41.730668 


E20 


47.203213 


44.789203 


43.111449 



Table 4: The lowest calculated energy levels for the cases k = 0, k = 1, and k = — 1 in domain walls 
dominated Universe (in all cases, A = —15). As mentioned in the text for every positive eigenvalue there exist 
a corresponding negative one with identical absolute value. 
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We found that, for all these cases the behavior of the expectation value of the scale factor is similar to ones 
depicted in Figs. [Tito 1131 and never tends to the singular point. 

5 Conclusions 

In this work we have investigated closed, flat, and open minisuperspace FRW quantum cosmological models 
(k = 1,0, —1) with perfect fluid for the radiation, dust, cosmic strings, and domain walls dominated Universes 
({a = 1/3, 0, —1/3, —2/3}, respectively). The use of Schutz's formalism for perfect fluids allowed us to obtain a 
Schrodinger-like WD equation in which the only remaining matter degree of freedom plays the role of time. We 
have used Spectral Method and obtained accurate results for the eigenfunctions and eigenvalues. Physically 
acceptable wave packets were constructed by appropriate linear combination of these eigenfunctions. The 
time evolution of the expectation value of the scale factor has been determined in the spirit of the many- 
worlds interpretation of quantum cosmology. Since the expectation values of the scale factors for the cases 
considered here never tend to the singular point, we have an initial indication that these models may not have 
singularities at the quantum level. The similar conclusions have been obtained on general grounds in |18j and 
for the radiation case in [22] . 
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